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■ Abstract 

We discuss the 2-point-particle-iiTeducible {2PPL) expansion, which sums bubble 
, graphs to aU orders, in the context of SU{N) Yang-Mills theory in the Landau gauge. 

' Using the method we investigate the possible existence of a gluon condensate of mass di- 

• ^ , mension two, ^AJJ A"), and the corresponding non-zero vacuum energy. This condensate 

^ ' gives rise to a dynamically generated mass for the gluon. 

M\ 

1 Introduction. 

Recently there has been growing evidence for the existence of a condensate of mass dimension 
two in SU{N) Yang-Mills theory with colours. An obvious candidate for such a condensate 
is I^A'^A'^. The phenomenological background of this type of condensate can be found in 

Q|2j|^. However, if one first considers simpler models such as massless theory or the 
Gross-Neveu model and the role played by their quartic interaction in the formation of a 
(local) composite quadratic condensate and the consequent dynamical mass generation for the 
originally massless fields 1^131121; it is clear that the possibility exists that the quartic gluon 
interaction gives rise to a quadratic composite operator condensate in Yang Mills theory and 
hence a dynamical mass for the gluons. The formation of such a dynamical mass is strongly 
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correlated to a lower value of the vacuum energy. In other words the causal perturbative 
Yang Mills vacuum is unstable. From this viewpoint, mass generation in connection with 
gluon pairing has already been discussed a long time ago in, for example, 0|Hlini|TO]. There 
the analogy with the BCS superconductor and its gap equation was examined. It was shown 
that the zero vacuum is tachyonic in nature and the gluons achieve a mass due to a non- 
trivial solution of the gap equation. Moreover, recent work using lattice regularized Yang 
Mills theory has indicated the existence of a non-zero condensate, (^^^A^^, [TT]. There the 
authors invoked the operator product expansion, (OPE), on the gluon propagator as well as 
on the effective coupling Og in the Landau gauge. Their work was based on the perception 
that, even in the relatively high energy region (~ lOGeV), a discrepancy existed between 
the expected perturbative behaviour and the lattice results. It was shown that, within the 
momentum range accessible to the OPE, that this discrepancy could be solved with a 
power correction*. They concluded that a non- vanishing dimension two condensate must 
exist. Further, the results of ^21 give some evidence that instantons might be the mechanism 
behind the low-momentum contribution to condensate. As has been argued in [51, only the 
low-momentum content of the squared vector potential is accessible with the OPE. Moreover, 
they argue that there are also short- distance non-perturbative contributions to (^^fij- 

It is no coincidence the Landau gauge is used for the search for a dimension two con- 
densate. Naively, the operator is not gauge- invariant. Although this does not prevent 

the condensate (AJ^) showing up in gauge variant quantities like the gluon propagator, we 



should instead consider the gauge-invariant operator {VT)~^ uimu J d'^x (^[f) ! where VT 
is the space-time volume and U is an arbitrary gauge transformation in order to assign some 
physical meaning to the operator. Clearly from its structure this operator is non-local and 
thus is difficult to handle. However, when we impose the Landau gauge, it reduces^ to the 
local operator yl^. Moreover, it has been shown that ^^^^ is (on-shell) BRST invariant 

|14[ I15 | [T6 l I17j. Another motivation for studying ^^^^ is the perceived connection between 
the gluon propagator and confinement. (See JH] and references therein.) More precisely, 
the gluon propagator exhibits an infrared suppression, as has been reported in many lattice 
simulations, jElEOll^ and using the Schwinger-Dyson approach, |221 EHl 01] • A dynamical 
gluon mass might serve as an indication for such a suppression. An attempt has already been 
made to explain confinement by a dual Ginzburg-Landau model or an effective string theory, 
in the Landau gauge, with the help of ^^^^ j25j . The fact that ^^^^ might be central to 
confinement, is supported by the observation that it undergoes a phase transition due to the 
monopole condensation in three dimensional compact QED [21 [Hj . 

From these various analyses the importance of ^^^^ must have become clear. Therefore, 
the aim of this article is to provide some analytical evidence that gluons do condense. To 
our knowledge, [22] is the only paper which effectively calculates ^^^^) without referring 

*The power correction due to the (G^^) condensate is too weak at such energies to be the cause of the 
discrepancy. 

^Although this equality is somewhat disturbed by Gribov ambiguities |13l . In this paper Gribov copies are 
neglected since we will work in the perturbative Landau gauge and sum a certain class of bubble diagrams in 
this particular gauge. It is a pleasant feature of the Landau gauge that (^^) can be given a gauge invariant 
meaning. In another gauge, the bubbles will no longer correspond to (^^) and the correspondence with 
( A?, ) is more of academic interest. 

\ 'V min 
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to lattice regularization. In j26j the standard way of calculating the effective potential for a 
particular quantity was followed, and all the problems concerning the fact that the considered 
quantity was a composite operator were elegantly solved. 

In a previous paper ^3 > we have discussed the 2PPI expansion for the Gross-Neveu model 
and found results close to the exact values for the Gross-Neveu mass gap and the vacuum 
energy. The 2PPI expansion does not rely on the effective action formalism of Instead 
it is directly based on the path integral and the topology of its Feynman diagrammatical 
expansion. In this paper we will discuss how to apply it to SU{N) Yang-Mills theories in the 
Landau gauge. Of course, it is not our aim to provide a complete picture of (^^fi^ but rather 
give further evidence for its existence since it lowers the vacuum energy. 



2 The 2PPI expansion. 

The SU{N) Yang Mills Lagrangian in d-dimensional Euclidean space time is given by 

= -^GI^^GI,^ + Cgauge+F.P. (2.1) 

where G^^, is the gluon field strength, 1 < a < N"^ — 1, Cgauge+F.p. implements the Landau 
gauge and its corresponding Faddeev-Popov part and ^ and ^ denote the ghosts and anti- 
ghost fields respectively. Issues concerning the counterterm part of (|2.1|) will be discussed 
later. First, we consider the diagrammatical expansion for the vacuum energy which we 
denote by E. As is well known, this is a series consisting of one particle irreducible, (IPI), 
diagrams. These IPI diagrams can be divided into two disjoint classes: 

• those diagrams not falling apart into two separate pieces when two lines meeting at the 
same point x are cut, which we call 2-point-particle-irreducible, {2PPI); (an example 
is given in Fig. 1) 

• those diagrams falling apart into two separate pieces when two lines meeting at the 
same point x are cut which we call 2-point-particle-reducible, {2PPR), while x is called 
the 2PPR insertion point; (an example is given in Fig. 2). 

We may now resum this series of 2PPR and 2PPI graphs, where the propagators are the 
usual massless ones, by retaining only the 2PPI graphs, whereby the 2PPR insertions, or 
bubbles, are resummed into an effective (mass) The bubble graph gluon 

polarization is then given by 

yiab _ _ 9_ 



+ 



feadfebc (^^A^) " (^^^"))] (2-2) 



where fate the SU{N) structure constants. We define the vacuum expectation value of 
as 

A = (a^^A-) . (2.3) 
The global SU{N) symmetry can then be used to show that 
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Figure 1: A 2PPI vacuum bubble. 




Figure 2: A 2PPR vacuum bubble, x is the 2PPR insertion point. 



Substitution of in ((T^ yields 

K = - S'^.'-^^'^X^^ (2-5) 
which results in an effective mass, m, running in the 2PPI propagators, given by 

-2 2 N d-l^ 

m = 9 ^^^A . (2.6) 

If we let E2PP1 be the sum of the 2PPI vacuum bubbles, calculated with the effective 2PPI 
propagator, then this E2PP1 is not equal to the vacuum energy E, because simply removing 
all 2PPR insertions is too naive. For instance, there is a double counting problem which is 
already visible in the 2PPR diagram of Fig. 2 where each bubble can be seen as a 2PPR 
insertion on the other one. However, we can resolve this ambiguity. A dimensional argument 
results in 

E = E2PP1 + cg^^^ (2.7) 

where c 7^ will accomodate the double counting. To determine the appropriate value of c, 
we use the path integral which gives 

1 

4- 



^-Ja^JadeU\AlA^^Al) . (2.8) 



The first two terms contribute unambiguously to the 2PPI part. For the last term, we rewrite 

' Ab AC Ad Ae\ _ / Ab ac\ / Ad Ae\ , / Ab Ad\ i ac Ae\ 



+ {A^Al) (AlAl) + (^AlAlA%A^)^ ^^^^ . (2.9) 
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Using (|2.4|) and the properties of the structure constants fate, we obtain 



dg^ Ag \\ ^ " " ^^ ''/2PPI 2g \ ^ '^/2PPI 

+ -fabcfade{AlAlAf^Al)^^^^ + 4 ^2 _ ^ ^ 
_ dE,ppj 1 N d-1 

From H2.7() . we derive 

dE_ ^ 9^2PPJ _^ 9^2PPJ dm^ ^ ct^ ^Icg^lS.— (2 11) 

dg^ drn} dg^ dg^ 

Combining (fTTUI) and (|TTT|) gives 

dE2PPi f N d-1^, 2 ^ d-ldA\ 1 TV d-1 



A + o^ - - 1_1 = 1 lA^ (2 12) 

drrf \N^-l d ^ N^-l d dgV 4 - I d ' ^ ' 

Then a simple diagrammatical argument gives 

dE2ppi A 



dfn'^ 2 



(2.13) 



which is a local gap equation, summing the bubble graphs into m^. Using this together with 
finally gives 



1 N d-1 
^ ~ ~ 4iV2_i ^ 

It is easy to show that the following equivalence holds 



(2.14) 



dE2PPi A dE , , 

= ^ ^ 7^ = . 2.15 
crfn 2 dfm^ 

To summarize, we have summed the bubble insertions into an effective (mass)^, m^. The 
vacuum energy is expressed by 

g2 AT d _ 1 

E = E2PP, - XiV^^^ ■ (2.16) 

We stress the fact that (|2.16p is only meaningful if the gap equation (|2.15p is satisfied. This 
means we cannot consider m or A as a real variable on which E depends. It is a quantity 
which has to obey its gap equation, otherwise the 2PPI expansion loses its validity. This also 
means that E{fn), or equivalently -E(A), is not a function depending on m (A), in contrast''' 
with the usual concept of an effective potential V{(p) which is a function of the constant field 
ip. 

In order to ensure the usefulness of the 2PPI formalism for actual calculations, we should 
prove it can be fully renormalized with the counterterms available from the original (bare) 
Lagrangian, ^2.1^ . However, it is sufficient to say that all our derived formulae remain valid 
and are finite when the counterterms are included. This also implies the 2PPI mass m is 

^V{(p) also makes sense if ^ 7^ 0. 
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renormalized and gives rise to a finite, physical mass^, ?7iphys- Furthermore, no new coun- 
terterms are needed to remove the vacuum energy divergences. The renormahzabihty of the 
2PPI expansion has been discussed in detail in 27 in the case of the Gross-Neveu model. 
Since the arguments for Yang Mills theory are completely analogous, we refer to j27j for 
technical details concerning the renormalization. 

Another point worth emphasising here, is the renormalization group equation, (RGE), 
for E. The first diagram of E2PP1 is given by the 0-bubble. Using the MS renormalization 
scheme in dimensional regularization in d dimensions, we find 
3 _ 1 / ^2 I d 

" = 4- 



167r2 



-m 



In 



2 

m 



1 



Ag^d-l N 



-m 



(2.17) 



Since E is a physical quantity, it should not depend on the subtraction scale /x. 
expressed by the RGE 



dE 
d^I 



d_ 

dp 



d 



2\ — 2 
TO 



d 



dm 



E 







where (3{g ) governs the scaling behaviour of the coupling constant 



dp 



(a 



09^ 



and k((7^) is the anomalous dimension of m? 




m 



705^ + 7i/ 



+ 



72/ 



+ 



The coefficients can be found in gHl Hi for /? and in [201 EHl EHl for 7^2 



This is 
(2.18) 

(2.19) 

(2.20) 
(2.21) 



/5o 



70 = ^ 



11 

y 

35 

6" 



N 
167r2 

N 
167r2 



(3i 



7i 



34 / iV 
Y U67r2 
449 / N 
H V 167r2 



/32 



2857 / N 



72 



When we combine all this information and determine p^ up to lowest order in g'', we find 



54 V 167r2 
75607 9C(3) 



864 



16 



N 
167r2 



(2.22) 



(2.23) 



tt2 



dp 







(2.24) 



Apparently, it seems that E does not obey its RGE. However, this is not a contradiction 
because of the demand that the gap equation ()2.15|) must be satisfied. The gap equation 
implies that In ^ oc ^ + constants. The consequence is that all leading logarithms contain 
terms of order unity. Hence, we cannot show that the RGE for E is obeyed order by order. 
The same phenomenon extends to higher orders. In other words, knowledge of p^ up to a 
certain order n, would require knowledge of all leading and subleading log terms to order n, to 
show explicitly that p^ = 0. We must therefore be careful not to interpret the non-vanishing 
of the RGE as a reason to introduce a "non-perturbative" /3-function, as is sometimes done, 



^"^phys ^® pole of the gluon propagator. 
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Figure 3: The first diagrams contributing to E2PP1 



3 Results. 

Up to 2-loop order in the 2PPI expansion (see Fig. 3), we find in the MS scheme 



E{A) 



3 g^N —2 27 g^N 



16 -I 



A 



fin 1^1^ 



64 167r2 iV2 
2 

+ 



1 



259 



In 





-i) 






3 /T 


1043 




32 



9 g^N^ A 
16 (167r2)2 _ 1 



891 



63 



+ S2 C(2) 

16 8 ^ 



where C('^) is the Riemann zeta function, 
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S2 



0(2 



0.2604341 



(3.1) 



(3.2) 



and Ci2{x) is the Clausen function. We have computed the relevant two loop vacuum bubble 
diagrams to the finite part using the massive gluon and massless ghost propagators which are 
respectively 



1 



p2 + 777,2 



and 



(3.3) 



in the Landau gauge. The expressions for the general massive and massless two loop bubble 
integrals were derived from the results of implemented in the symbolic manipulation 



language Form, ISB]- It is easy to check that solving the gap equation |^ = 0; with fj? set 

7;2/ 



equal to | ^i^i ^ to kill potentially large logarithms, gives no solution at 1-loop or 2-loop, 

apart from the trivial one A = 0. This does not imply (^^^^ does not exist but that the 

MS scheme might not be the best choice for the 2PPI expansion. To address this we first 
remove the freedom existing in how the mass parameter A is renormalized by replacing A by 
a renormalization scheme and scale independent quantity A. This can be accommodated by^ 



A 



/(/) 



A 



with 



-Of 



(3.4) 



(3.5) 



^Barred quantities refer to the MS scheme, otherwise any other massless renormahzation scheme is meant. 
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A change in massless renormalization scheme corresponds to relations of the form 



A 



bog' + hg" + ■ 
A(l + dog"^ + dig^ + • 
/(5')(l + /o5' + /i/ + ---) . 



(3.6) 
(3.7) 
(3.8) 



With these, it is easily checked that (|3.4)) is renormalization scheme and scale independent. 
The explicit solution of 1)3. 5|1 reads 



7{f) 




Pi 



/3i7i 72 

Pi Po 



+ 0{g 




Pilo 7i^ ^ 
Po. 



P} 







(3.9) 



Since the gap equation is still a series expansion in and we hope to find (at least qualita- 
tively) acceptable results, jg;^ should be small. We will therefore choose to renormalize the 
coupling constant in such a scheme so that E is of the form 



E 



3 
16 



iV2 - 1 



167r2 1 



{I 



N 



E^L + E^L'^) + ...\ (3.10) 



where 



L 



In 



Ha 



A 



(3.11) 



Otherwise, we remove all the terms of the form (^f^^ x constant, and only keep the terms 

that contain a power of the logarithm L. This is always possible by calculating the MS 
value of E as in 1)3. 1|) and using 1)3.6(1 to change the coupling constant renormalization by 
a suitable choice for the coefficients bi. In other words the 1-loop MS contribution to E 
allows one to determine bo- Once bo is fixed, the 2-loop MS contribution to E can be used 
to fix bi, and so on for the higher order contributions. We note that the gap equation 
is translated into ^ = since oc A oc A. In this gap equation, we will set 

3 f„2\^ 2l3o 



4 vy y " so that all logarithms vanish. In other words L = 0. Notice that one 



cannot set /i 



2/3o 



N 



4 \y J " iv^irT expression ()3.1U|) and use the RGE for E to sum 

the logarithms when ^ = is solved. As already explained in the previous section, the RGE 

dA ^ 

for E is not obeyed order by order, see also [^. Once a solution A^, of the gap equation is 
found, then we will always have 



_3_ 
16 



70 



N 

iV2 - 1 



A 



(3.12) 
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If the constructed value for is small enough, then we can trust, at least qualitatively, the 
results we will obtain. Now we are ready to rewrite H3.1() in terms of A. After a little algebra, 
one finds 



16 



with 



iV2 - 1 



-1 + 



3N' 
2567r4 

9iV 



9N / 5 \ 

/ 9N 



259 ^ 31 ^ 
C3 + —L - —L^ ) - 26oCi - cf - 2c2 + 



C5 + 



647r2 



6 



- + L6o + ci+2-- + Lci + 6o l--^ 



6 



V647r2 I g 



70 



- + L - 2ci C4 



2/3o 



9' + 0(d 



ci 

C2 
C5 



1 / /?i7o 7i 



2 V f^l /3o 



1043 63 



V 

, 891 
32 8 ^(') + T^^^ 



^}_f Pill 72 



4 V Pi Po 



u2 70 



^ _ To 



2/3o V Po 



Next, we determine bo and bi so that H3.13() reduces to 



(3.13) 

(3.14) 

(3.15) 

(3.16) 
(3.17) 
(3.18) 



E 



i_2a 



Po ^2 



N 



m - 1 



9A^ 

+ ^(^o + ^4 + 2ci)L 



2 9iV , 4 

-^ + ^64^^ + ^ 



+ 5 



3iV^ /259^ 31^. 



2567r4 V 



(3.19) 



We find that feo is 



409 N 



We do not list the value for bi since it is no longer required. From the /3-function we find the 
two loop expression for the coupling constant is 



/3i In In 



/3oln^ Poplln^ii^ 



(3.21) 



where A is the scale parameter of the corresponding massless renormalization scheme. We 
will express everything in terms of the MS scale parameter Aj^. In jSl], it was shown that 



to 

A — A P 2/3o 



(3.22) 
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We will also derive a value for the (^G^j/) condensate from the trace anomaly 



G 



m 

2<? 



(3.23) 



This anomaly allows us to deduce for = 3 the following relation between the vacuum energy 
and the gluon condensate 

32 

— ^vac • (3.24) 



At 1-loop order, the results for = 3 are 
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167r2 



a. 



G 



TT 



1-loop 



1-loop 
1 1-loop 



1-loop 



1.004A^^ 
-0.0074A|j^ 
0.02Aij^ - 



233MeV 
« - 0.00002GeV'^ 
0.00007GeV^ 



(3.25) 

(3.26) 
(3.27) 
(3.28) 



while the scale parameter jf' (184MeV)^. We have used Aj^ ^ 233MeV which was the 
value reported in We see that the 1-loop expansion parameter is quite large and we 

conclude that we should go to the next order where the situation is improved. We find 



167r2 



2-loop 



2-loop 
vac 1 2-loop 



0.131 



2.3Aj^ « 536MeV 



-0.63Aijg ^ - 0.002GeV^ 



a 



IT 



2-loop 



1.84A|^ 



0.005GeV^ 



(3.29) 

(3.30) 
(3.31) 
(3.32) 



with ]f (347MeV)2. Although there is a sizeable difference between 1-loop and 2-loop 
results, the relative smallness of the 2-loop expansion parameter, indicates that the 2-loop 
values are qualitatively trustworthy. It is well known that in order to find reliable perturbative 
results, one must go beyond 1-loop, and even beyond 2-loop approximations. Therefore, one 
should not attach a firm quantitative meaning to the numerical values. Let us compare 
our results wi th wh at was found elsewhere with different methods. A combined lattice fit 

resulted in (^m) ~ 1.64GeV, We cannot really compare this with our result for 

since the lattice value was obtained with the OPE at a scale fx = lOGeV in a specific 
renormalization scheme (MOM). However, it is satisfactory that (|3.30jl is at least of the same 
order of magnitude. More interesting is the comparison with what was found in 26_ with the 
local composite operator formalism for (^^fij- In the MS scheme at 2-loop order, it was found 

that — T 



0.141247 while ^vac - 
results. An estimate of the tree level gluon mass of 



0.789A|-j^ which is in quite good agreement with our 



SOOMeV was also given in which 
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compares well with the lattice value of ~ 600MeV of [351 136j . With the 2PPI method, one 
does not really have the concept of a tree level mass. Instead, the determination of Wpjiyg 
would need the calculation of the highly non-trivial 2-loop 2PPI mass renormalization graphs 
which is beyond the scope of this article. 

In conclusion we note that the perturbative Yang-Mills vacuum is unstable and lowers 
its value through a non-perturbative mass dimension two gluon condensate (^^fij- We have 
omitted quark contributions in our analysis but it is straighforward to extend the 2PPI 
expansion to QCD with quarks included. Indeed an idea of the effect they have could be 
gained by an extension of PH] . 
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